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We develop a general theory of non-equilibrium states based on the Keldysh formalism, 
in particular, for charged-particle systems under static uniform electromagnetic fields. The 
Dyson equation for the uniform stationary state is rewritten in a compact gauge-invariant 
form by using the Moyal product in the phase space of energy-momentum variables, whcich 
originally do not commute in the case of the conventional operator algebra. Expanding 
the Dyson equation in electromagnetic fields, a systematic method for the order-by-order 
calculation of linear and non-linear responses from the zeroth-order Green's functions is 
obtained. In particular, we find that for impurity problems, up to linear order in the electric 
field, the present approach provides a diagrammatic method for the Stfeda formula. This 
approach also generalizes the semi-classical Boltzmann transport theory to fully quantum- 
mechanical and/or multi-component systems. In multi-component systems and/or for Hall 
transport phenomena, however, this quantum Boltzmann transport theory, constructed from 
the anti-symmetric combination of two different representations for the Dyson equation, does 
not uniquely specify the non-equilibrium state, but the symmetric combination is required. 
We demonstate the formalism to calculate longitudinal and Hall electric conductivities in an 
isotropic single-band electron system in the clean limit. It is found that the results are fully 
consistent with those obtained by Mott and Ziman in terms of the semi-classical Boltzmann 
transport theory. 

§1. Introduction 

Theory of non-equilibrium states is one of the most fundamental and important 
topics in statistical and condensed-matter physics. For instance, applying an electric 
field drives a charged system off the equilibrium. Electron transport properties in 
a condensed matter under an applied electric/magnetic field are conventionally de- 
scribed in terms of the semi-classical Boltzmann transport theory EJS Rigorous for- 
mulations based on quantum mechanics have also been developed by several authors. 
The theory of the linear response of equilibrium systems to external perturbations 
was established by KubcP and have also been formulated in various manners) 4 ' 
In this case, the linear responses are described by correlation functions of the equi- 
librium state. To go beyond the linear response in the Kubo formalism, one needs to 
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expand the correlation function further in the external perturbation. From this view- 
point, Fukuyama and coworkers performed a diagrammatic calculation of the Hall 
effect in metals in the presence of impurities which has also been generalized for 
Fermi liquids P^J 

The KeldyshpJ or Baym-Kadanoffp^ formalism is a powerful theoretical frame- 
work to treat non-linear and/or non-equilibirum phenomencP^'^'^ by directly 
handling the two-point Green's functions of all the independent analyticities. In 
particular, in the Wigner representation in terms of the center-of-mass time-space 
coordinate and the relative canonical energy-momentum p^, each product ap- 
pearing in the Dyson equation for a uniform equilibrium, i.e., X^-independent, sys- 
tem is replaced by the so-called Moyal produdP^ in this "phase" space, in order to 
fully describe the convolutions of X^-dependent Green's functions.^ 

It is known that taking an antisymmetric combination of two different rep- 
resentations for the Dyson equation and then performing the gradient expansion, 
the Keldysh formalism gives a fully quantum-mechanical generalization of the semi- 
classical Boltzmann transport theory,^ EJEJ which is often referred to as the quan- 
tum Boltzmann transport theory. In particular, for uniform stationary states under 
constant, i.e., uniform and static, electromagnetic fields, the Green's functions can be 
expressed as apparently gauge-invariant functionals of the energy-momentum vari- 
able 7r^ = p 11 — qA^(X). Here, q is the electric charge of the particle and A^(X) is 
the scalar and vector potential of the electromagnetic field. This ir^ is the Wigner 
representation of the mechanical or kinetic momentum^) or the covariant derivative, 
which should not be confused with the canonical momentum p^ . 

For open systems, the Keldysh formalism is combined with the Landauer-Buttiker 
formalism^ by imposing a certain boundary condition, or one can deal with 
the time-dependent Schrodinger equation to observe the time evolution of a non- 
equilibrium stated In closed systems, however, it is necessary to develop system- 
atic ways of solving non-linear and/or non-equilibrium problems without any semi- 
claasical approximation. Not even the linear response in multi-component systems 
and Hall effects can be described solely by the quantum Boltzmann equation, even 
though it is correctly generalized for these cases. Instead, in addition to the quan- 
tum Boltzmann transport equation, one must directly handle the Dyson equation 
under the external fields in multi-component systems™ This point was recently 
demonstrated in the context of the spin Hall effect 211 ' 22 ' and the anomalous Hall ef- 
fect It is remarkable that the semi-classical Boltzmann equation is not sufficient 
for explaining the Hall transport phenomena. 

In this paper, we develop a generic theoretical framework for systematically 
calculating the linear and the non-linear responses to an applied uniform and steady 
electromagnetic field. We obtain a compact form for the Dyson equation in terms of 
the Moyal product in the tt space in the presence of the electromagnetic field tensor 
F^ u , which re-expresses the non-commutativity of the components of the mechanical 
momentum in the conventional operator algebra. This yields a master equation 
for solving non-linear and/or non-equilibrium but uniform stationary problems. Up 
to the linear response to the electric field, we obtain an expression for the electric 
conductivity in the form of the Stfeda formulaP^ which was originally intended 
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for non-interacting electrons. Now we reveal the diagrammatic technique in the 
Streda formula for more generic cases. Generalization to non- uniform and/or time- 
dependent problems is also straightforward. 

We demonstrate the formalism by applying it to the conductivity tensor, in- 
cluding the longitudinal conductivity and the transverse Hall conductivity in an 
isotropic single-band electron model in the presence of impurities. In the clean limit, 
our results coincide with those of the semi-classical Boltzmann transport theoryP^ 
We also found the correspondence to diagrammatic calculations of vertex correc- 
tions due to the impurity potential in terms of the Kubo formula in the case of the 
linear-response longitudinal conductivity.^ 

Major advantages of the present formalism are the following. 

1. The Green's functions can be described in terms of gauge- invariant functionals 
of the energy-momentum 7r M = — qA^{X) in an arbitrary order in the elec- 
tromagnetic fields. Therefore, the gauge invariance of physical quantities after 
integration over n is evident from the beginning. 

2. It allows for the electric current to be expressed as a sum of the equilibrium 
thermodynamic part, which is proportional to the Bose or Fermi distribution 
function /=p(e), and the nonequilibrium parts given by the energy (e) derivatives 
9™/zp(e) (m = 1, • • • , n) in the case of the n-th-order response. This enables 
us to carry out a diagrammatic calculation of the Streda formulc!^ and its 
generalization to non- linear responses. 

3. It is easy to treat multi-band/component systems. In general, not only the mag- 
netic field but also the electric field alters the retarded and advanced Green's 
functions, even in the absence of an interaction. This is beyond the scope 
of the conventional quantum Boltzmann theory, but it can be systematically 
calculated in the present approach. 

This paper is organized as follows. In £J2J we briefly introduce the Keldysh formal- 
ism for non-equilibrium Green's functions. Then, we formulate a generic framework 
to systematically calculate the bulk quantum distribution function in response to 
the uniform static electromagnetic fields for generic multi-component bosonic and 
fermionic systems. Sections are devoted to explicit applications of this generic 
formalism. In the linear response theory is derived. In it is generalized to the 
cross terms of the electric and magnetic fields. Then, in we study a simple model 
of an isotropic single-band electron system as a concrete application of the present 
formalism. The equivalence to the Kubo formula and the Streda formula and the 
merits of this approach in practical calculations are demonstrated. We have also 
obtained results for the longitudinal and the Hall conductivities that are consistent 
with the semi-classical ones obtained by MottP and ZimanP Section H3 is devoted 
to discussion and conclusions. 

§2. Keldysh formalism for non-equilibrium Green's functions 

To consider responses of a charged-particle system to external electromagnetic 
fields, we employ the Keldysh formalism^ for non-equilibrium Green's functions. 
Here, in addition to the retarded and advanced Green's functions, one introduces the 
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Keldysh or the lesser Green's function. In particular, the lesser Green's function plays 
the role of the quantum distribution function. These non-equilibrium and/or non- 
linear Green's functions can be calculated at least by performing an order-by-order 
expansion in the electromagnetic fields. 

2.1. Model 

We consider a generic multi-component model of bosons or fermions with an 
electric charge q whose Hamiltonian is given by 

H t ot (x,p) = H (p) + V (x), (2-1) 

with the space coordinate x = (x,y,z) and the conjugate momentum p = —ih'V x . 
Here, Hq(p) and Vq(x) denote the non-interacting and the impurity-potential parts 
of the Hamiltonian, respectively. Henceforth, a variable with a hat represents a 
matrix in the space of internal degrees of freedom including spins. 

Next, we introduce the external electromagnetic field F^ u = d XfJ/ A u (x)— d Xv A^{x). 
Henceforth, x M = (t, x) and x^ = (— i, x) denote the time and space coordinates, 
d xf j, = (d t ,V x ) and d x = (— d t , V 2 ) the derivatives, and A^{x) = (<p(x),A(x)) and 
Afi{x) = (—(p(x), A(x)) the sets of scalar and vector potentials yielding the electro- 
magnetic field tensor F fMU (x) = d x nA u (x) — d x ^A^(x). Then, Eq. (|2Tj) is modified by 
A(x) as 

H tot (x,p) = Hl ot {x,p- qA{x)) + q<t*(x) 

= H {p-qA{x)) + %{x)+qcj>{x). (2-2) 

Although we have included neither the Coulomb interaction among the electrons 
nor the phonons into the Hamiltonian, the following formalism based on the Dyson 
equation holds even in the presence of the Coulomb interaction and the phonons 
if the corresponding self-energy corrections are properly included according to the 
Feymann ruleP^* 

2.2. Green's functions 

The retarded, advanced, Keldysh and lesser Green's functions are defined as 

^xi;x 2 ) = -i6(ti-t 2 )(\Va 1 (xi)i 1 Pi 2 (x2)\_)i (2-3a) 



\a 2 ) 
;a 2 ) 



I > (ai 

{g k \ 

y > (<»l 

G<\ r (x 1 ;x 2 ) = ^ihl 2 (x 2 )^ ai (.x 1 ) 

J (ai;o 2 ) \ 



(xr,x 2 ) = i0(t 2 -*i)([lP r ai(iCi),*i,(a?2)] ), (2-3b) 
(xi;x 2 ) = -i(\^ ai (xi),^a 2 (x 2 )}), (2-3c) 



= l - (-G R (x l]X2 ) + G A ( Xl ;x 2 ) + G K (x i; x 2 )) , (2-3d) 

respectively. Here, the upper and lower signs correspond to bosons and fermions, 
respectively. Wa m {x m ) and ^ Qm (x m ) represent the creation and annihilation opera- 
tors of a boson or a fermion with the internal degree of freedom labeled by a m at 
the time-space coordinate x m . 
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In the present approach, the lesser Green's function represents the electron dis- 
tribution function (multiplied by 2iri) even in the presence of the external fields and 
is defined for the full time-space variables. This is in sharp contrast to the situa- 
tion in the semi-classical Boltzmann transport theory, where only its simultaneous 
component with t\ = t 2 is considered. In the semi-classical case, quasiparticles are 
implicitly assumed to be well-defined so that the broadening of the peak in their 
spectra is negligibly small compared with the other energy scales. By contrast, the 
present approach is applicable even in the case of a strong scattering amplitude for 
which the usual elastic approximation and thus the semi-classical Boltzmann trans- 
port theory become invalid. 

2.3. Dyson equations 

In the conventional Larkin-Ovchinnikov representation, one introduces the fol- 
lowing matrix-form Green's function and self-energy: 

G R {x 1 -x 2 ) G K (x 1 ;x 2 
G A {xi;x 2 ) 

S R (xi;x 2 ) U K (xi;x 2 
Z A (x 1 ;x 2 



G{ Xl -x 2 )=[ " ^' X2) Xl " . , )> ( 24 ) 



E(x 1 ;x 2 )=( " ^ X2 > % A )y_-: 2 ( ). (2-5) 



Here, 8(xi — x 2 ) represents 5(t\ — t 2 )5(xi — x 2 ), and S R ^ A \xi;x 2 ) and £ K ( < >(x\;x 2 ) 
denote the retarded (advanced) and Keldysh (lesser) self-energies, respectively. In 
the present paper, however, we further perform an orthonormal transformation in 
terms of 



as 



/■t, , r\ru i/o- 1 i G (^l!^) 2G < (xi;x 2 ) i ,.-)-\ 

G(x 1 ;x 2 ) <- Q G{xi;x 2 )0 = \ - A / / , (2-7) 



G A ( Xl ;x 2 

fv \ /~« w \r>-l ( ^ R (xi)X 2 ) 2E < (xi;x 2 ) \ , , 

E{ Xl ,x 2 )^OE{x x] x 2 )0 =^ Q E A {xi . X2) )■ (2-8) 

This representation has the advantage of allowing the direct calculation of the lesser 
component without any manipulation of the Keldysh components. 

According to Eq. (|2-2|) . we can write the equations of motion or the Dyson 

equations for the Green's function in terms of the gauge covariant derivative D Xm/i = 

dx m)1 ~ iqA^(x m )/h and its conjugate D* m ^ = d Xm ^ + iqA^(x m )/h as 

iHD tl G(xi;x 2 ) - Ho(-ihD Xl )G(xi;x 2 ) - (£_* G)(xi;x 2 ) = 5(xi - x 2 ), (2-9a) 
-i6(x 1 ;x2)hDt a -G( Xl ;x 2 )H (ihD* X2 ) - (6 * E)( Xl ; x 2 ) = 8{x x -x 2 ). (2-9b) 

Henceforth, d and d denote the derivatives operating on the right-hand and left- 
hand sides, respectively, and * represents the convolution, defined by 

(Z*K)(xr,x 2 ) = j dx 3 f(xi;x 3 )S(x 3 ;x 2 ), (2-10) 
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where J dx represents fdtf dx. The Dyson equations (|2-9a[) and (|2-9b|) can also be 
written in the forms 



((G^- 1 -r)*G)(x 1 ;x 2 ) = 5(x 1 
(G*(G(°)- 1 -^))(x 1 ;x 2 ) = 5(x 1 



where 



G(°)- 1 (x 1 ;x 2 ) 



ihD tl - H (-ihD xl 



-ihDl 



H (ihD* ) 



5(xi 
5(x\ 



X2) 
X2) 



(2-lla) 
(2-llb) 

(2-12a) 
(2-12b) 



It is known that the self-energy 17 can be constructed from a Feynman rule^J similar 
to that in the conventional Green's function method for G R ' A and E R,A . 

2.4. Wigner representation and non-commutativity as the Moyal product 

We change the set of variables (2152:2) to (X;x) with = (T,X) and x^ = 
(t, x) being the center-of-mass and relative coordinates, respectively: 

X = ( Xl +x 2 )/2, (2-13a) 
x = x\ — X2- (2-13b) 



We proceed to the Wigner representation in terms of (X;p) by means of the 
Fourier transformation from x^ to p^ = {p\^ — p2 IM )/2 = (e,p) or from x^ to p^ = 
(Pin ~ P2 M )/2 = (—e,p). It is well known that the convolution given in Eq. (|2-10|) 
can be rewritten in a compact form as a function of X^ and p^; 



(±*i_)(X;p) =f{X-p)e^ a ^ d ^- d ^ dx ^S{X-p), (2-14) 

which is explicitly proven in Appendix EJ Henceforth, contraction of indices for 
four-component vectors should be understood. This takes the form of the so-called 
Moyal product or the star product of T(X;p) and S(X;p) in the space of X^ and 
p^, which fully represents the non-commutativity [xf,p^]- = ih. Recently, this has 
also been utilized in the field of high-energy physics in the context of the D-brane 
in string theorypZJ and the associated non-commutative quantum mechanics^ has 
been intensively investigated. 

Now, in order to express the Dyson equations in the Wigner representation, 
we have only to calculate &°^~ 1 (X;p). For this purpose, it is useful to define the 
quantity associated with the covariant derivative 



-ihD Xl ii8{x\ 



x 2 ) = ihD* ^5{x\ - x 2 



(2-15) 



so that GW-^xi; x 2 ) can be represented in the form of the convolutions of tt^ (x\ \x 2 ). 
The Wigner representation of Eq. (|2-15|) associated with the covariant derivative is 
obtained as 



7r fl (X;p)= fdxe-^' h K^ 



x + x -,x- x - 

2' 2 



Vtl - qA^X) 



(2-16) 
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using integration by parts. This new variable TT^(X;p) defines the mechanical or 
kinetic energy-momentum in the Wigner representation, and it is clearly distinct 
from the canonical one, p^. 

Finally, the Wigner represetation of G^ ) -1 can be obtained in a compact form 
by using the Moyal * product as 



G^- 1 (X;p)=G ( S- 1 (7r(X;p)), 



where 



4° ) - 1 ( 7 r)=7r - J ff Q *(7r) 



and Hq^(tt) in the presence of electromagnetic fields are generalizations of 



G, 



(0)-l 



(p) 



H (p) 



(2-17) 
(2-18) 

(2-19) 



and Hq{p) defined in the absence of the fields. They are obtained by replacing all the 

products in G^ 1 (p) and Ho(p) with the Moyal * products defined by Eq. (|2-14l) 
and also p with tt = p — qA(X). Therefore, in the Wigner representation, the Dyson 
equations ()2Tla|) and (|2-llb|) are rewritten as 



7T°(X;p) - H *(tt(X; P )) - E(X;p) *G(X;p)=l, 



G(X;p) * ir°(X;p) - H ,(tv(X; P )) - E(X;p) 



1. 



(2-20a) 
(2-20b) 



2.5. Gauge invariance in the case of constant electromagnetic fields 

Here, it is convenient to decompose the electromagnetic potential into two parts; 
A(X), which just contributes to uniform static electromagnetic fields, = d x nA v (x) 
— d x v A^{x) = constant, and the other non-uniform and/or dynamical part, 5A(X) = 
A{X) — A(X). Similarly, we introduce 



^(X;p)=p»-qA(X). 



(2-21) 



Now, we change the set of variables from (X;p) to (X\tt). Then, the Moyal product 
given by Eq. (|2T4I) is transformed into the partly gauge-invariant form 

(f*S)(X;p) = (f>H)(X|7f) = f(X\n)e%^ x ^^-^^ x » +q ^^*" ) S{X\n). 

(2-22) 

In this (^|vr) representation, the Dyson equations take the forms 



7T°-q&f>(X) - H *(n - q5A(X)) - E{X\tx) *G(X\tx) = 1, (2-23a) 



G(X\tt)* tt° -q&j)(X) - H *(n - q 5A(X)) - S(X\tt) 



(2-23b) 



In the rest of this paper, we restrict ourselves to the case of a uniform static 
applied electromagnetic field F^ u = d x ^A u — d Xv A tl = constant, i.e., A{X) = A(X), 
or equivalently, 5A(x) = 0. Now we do not have to distinguish between tt and tt. In 
the case of the constant electromagnetic fields, the Moyal * product in the (^|vr) 
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space, which is defined by Eq. (J2-22J), gives the gauge-invariant representation of the 
conventional Moyal * product given by Eq. (|2-14() . Namely, like Y(X\ir) and 5(X\ir), 
their convolution has only an implicit dependence on the choice of the gauge through 
7r^ = — eA^i^X), and hence the gauge invariance after the p or tt integration is 
evident. This is consistent with Schwinger's idea to extract gauge-invariant results 
by using gauge-covariant quantities J2DJSDE3 since TT^{X]p) is constructed from the 
gauge-covariant derivative via Eqs. (|2-15|) and (|2-16|) . Furthermore, in the present 
case of constant electromagnetic fields, imposing the fully symmetric combination of 
the Moyal products among tt for the Hamiltonian Hq^(7t), as required for quantum- 
mechanical problems, the * products implicitly appearing in Hq(tt) are reduced 
to conventional products, owing to the antisymmetric property F^ u + F V)1 = 0. 
Therefore, we arrive at the following gauge-invariant forms of the Dyson equations 
under a constant electromangeitc field: 

tt -H (ir) -E{X\tt)\ *&(X\ir) = l, (2-24a) 
G(X\ir) * \ir° - H (tv) - £( x k)] = 1 - (2-24b) 

Furthermore, as long as we restrict ourselves to a uniform, stationary solution, 
Eq. (|2~2l?|) reduces to 

(Z>S)(t) = f{TT)e^ F ^^« v S{TT). (2-25) 

This, again, takes the form of the Moyal product, but now it translates the non- 
commutativity among the components of the mechanical energy-momentum opra- 
tor tt^ = p^ — eA IJ '(xi), which is introduced by the physical electromagnetic field 
hqF^ ^ 0. Then, using this Moyal product defined by Eq. 1)2-25(1 . the Dyson equa- 
tions (I2^4all and (!2-24bl) can be expressed in the gauge-invariant forms 

(tt -Hq(tt) -|Z(7r))*6(7r) = 1, (2-26a) 
G(tt) * (tt° - H {-k) -JZ(tt)) = 1. (2-26b) 

These two Dyson equations guarantee that for a given 27, the solution for G is unique. 
Therefore, one can define the inverse with respect to the Moyal product as 

G(tt) = (tt - H (tv) - J2(tt)) ~* . (2-27) 

2.6. Self-energy due to the scattering by dilute impurities 

Here, we consider the self-energy arising only from the scattering by dilute im- 
purities. For a general impurity potential of the form 

V( x ) = ^2,u(x-Xi), (2-28) 

i 



the simplest self-consistent Born approximation reads 

dir 

(2ttH) 



27(tt , tt) = n imp / 7 -^_u 7T _ 7T ,G(7T°, Tr')^-™ (2-29) 



J dxe~ iq - x u{x). (2-30) 
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When the potential is given by a (^-function, i.e., 

u(x) = u5(x), 

U q = U, 



(2-31a) 
(2-31b) 



it is easy to construct the self-consistent equations in terms of the T-matrix approxi- 
mation, taking account of an inifinite series consisting of the first and all the multiple 
Born scattering amplitudes for the self-energy as 



JZ(tt°) = n imp f (vr°), 
T(7r°) = u(l -g{TT°)u 

d7V 



(2-Khy 



G(VT U ,7T). 



(2-32) 
(2-33) 

(2-34) 



2.7. Expansion in the electromagnetic fields 

In most cases, except for the simplest problems, like the Zener problem,^ the 
Dyson equations l)2-26a|) and ()2-26b|) cannot be directly solved in the context of 
non-commutative quantum mechanics using the Moyal product. Therefore, the per- 
turbative treatment is an important and powerful tool. 

Dyson equations in the form of Eqs. (|2-26j) can be expanded in hqF^ u as 



7T° -H (ir) -£(tt)) G(tt) - 1 



Ei(''-».w-a-)) n 



n=l 



ihqF^ 



d d n "i G(tt) 



\i=l 



(2-35a) 



G(vr) (vr° - H (tv) - £(*-)) - 1 

n=l ' \i=l / 



n=l 

Finally, we obtain 

G(tt) = 4(7r) [l + ( r(vr) - r (vr) ) G(tt) 
1 



(2-35b) 



[ [ 7) d ^ d G(tt) 



71=1 



vi=l 



(2-36a) 



[i + G(tt) (M*)-Za{< 



71=1 



'9^5^ vr -F (7r)-r« 



Vi=l 



G (vr). 



(2-36b) 
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Henceforth, quantities with the subscript "0" denote those defined in the absence of 
external electromagnetic fields. This notation is introduced in order to distinguish 
them from those in the presence of the fields. 

We note that a non-linear generalization of the quantum Boltzmann equation in 
the steady, uniform case can be derived from Eqs. (|2-35f) as 



oo 1 

T- 

n=l 



ihqF^ 



\i=l 



7T° - Hn(-n) - ZYvr' 



n 



(2-37) 



in the Keldysh space. It is clear that this quantum Boltzmann equation cannot 
determine the retarded and advanced components G R ' A (tt) and E R)A (ir), as readily 
verified for single-component, i.e., scalar, cases. Furthermore, in multi-component 
systems, it does not uniquely determine G < (ir) and £ < (tt) either, as we will show 
later. 

Equations (|2-36b|) and (|2-36b|) have the solutions of the forms 



oo -. / n 

n=l ' \i=l 

oo 1 / n 

n=l ' \i=l 



2 



y 

— p,\Vl,— ,/lnfi 



(vr). 



(2-38) 
(2-39) 



We can calculate G(n) and by equating the coefficients of each order in F^ v in 
Eq. (|2-36b|) or (|2-36bj) and, for instance, imposing that G^,... ^ nVn and K^u u - ,p, nVn 
are symmetric under the permutation {fiiVi) <-> {^jVj) and anti-symmetric under the 
permutation //j «-»• Vi. 

In particular, in the case of impurity problems, we obtain 



y 

id-mux,— ,fl; 



dir' 



(2nh) 



(2-40) 



in the self-consistent Born approximation for the general impurity pontential given 
in Eq. ffigty . and 



( 7 r°) = n ini pio(- )EE E' II 

s 1=1 mi,— ,m.i_i j=l 
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(mj — Wlj-l)! -/ 1 s(m j _ 1 + l) l 's(m j _ 1 +l) ! ->/'s(m i } 1 's(m : i) 



(vr ) toiA 



9 (tt°) 



d7T 



(2vr^) 



d^-Ml "1,- ^n". 



(7r°,ir) 



(2-41) 
(2-42) 
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for the T-matrix approximation with the <5-function impurity potential given in 
Eq. (|2 31f) . where m(0) = and m(i) = n, and Yj m (i),— m(i-l) anc ^ denote 
the summations over m(l), • • • , m(i — 1) under the condition m(j) < m(j') for j < j' 
and over all the permutations among 1, • • • , n, respectively. 

2.8. Electric current 

Here, we summarize the formulation of the electric currents in the presence of 
constant electromagnetic fields. 

The component-dependent particle number and current density at a given time- 
space coordinate X 1 *, pP{X) = (n(X), j (X)), are expressed in terms of the quantum 
distribution function G < as 



n{X) 



-iCT 



X + -;X 
2' 



x->0 



°° dTT° 

.™ 27ri 



0l d 7V 

(2irh) d 



d d 7T 1 



_oo27ri J {2ith) d 2 



*o(7r),G < (X| 7 r) 



G < (X\tt), (2-43) 



(2-44) 



with the velocity 



v(p) 



1 

ih 



x,H tot (x,p) = V p H(x,p) = V n H (7v) = v (tt) 



(2-45) 



which satisfies the desired condition of the equation of continuity for the total particle 
number, 

d x ^{X) = 0, (2-46) 

with the total number and current density p^{X) = Tr/5 M (X). 

Finally, the uniform, static electric currents are given by the average over both 
time and space 



J = q lim lim 

H T^oo V^oo TV 



T/2 

i 

-T/2 



dT' / dX'Tr 



v 



(2-47) 



It is easy to see that Gq(tt) never contributes to the current J without breaking the 
inversion symmetry in the absence of an electric field. 

2.9. Zeroth-order equilibrium properties 

Lastly, for later use, we summarize the well-known equilibrium properties in the 
absence of external electromagnetic fields: 

-i 



e-H (p)-S^ A \p) 



Z£{p) - Z*(p)) f T (e), 



(2-48) 
(2-49) 
(2-50) 



with the Bose and the Fermi distribution functions /-(e) = (1 — coth(/3(e — /u)/2))/2 
and f+(e) = (1 — tanh(/3(e — / u)/2))/2, respectively. 
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The self-energy is obtained as 



^ (A) ie) = n 



imp 



dp 



(2irhy 



(2-51) 



in the Born approximation for the impurity potential given by Eq. Q2-3UJI . and as 

Z^ A \e) = n imp fK {A \e), (2-52) 



u 



(2irk) 



(2-53) 
(2-54) 



in the T-matrix approximation for the <5-function potential given by Eq. ()2-31(l . 

§3. Linear response theory 

We proceed to the linear response of the Green's functions to the external fields 
F^y. In particular, we must calculate G< u to consider the linear response of various 
quantities. For this purpose, we linearize Eq. ()2-36b|) or (|2-36b|) in F^ u to obtain 

<^»=G (7T) 



Go(vr), 
(3-1) 



or 



G E (ir) = Go(tt) 



&W+2^ (t))Go(t)(V.Go (tt)) 

- (V^V^W^o^V)))! GoW. (3-2) 



^b(vt) + ^(V.G '(tt)) x G (vr)(V.G '(tt)) 



G B (7r) = G (vr 

Here, we have used the identity 

d n nG (ir) = -G (ir) (d^d 1 (tv) \ G (ir), 

or equivalently 



G (7r). (3-3) 



(3-4) 



^R(A) 



iR(A)-l, 



(3-5a) 



d^G^i-rr) 



- G r (tt) (flU^W) G£(tt)] / T (vr°) 
+ (^(vr)-^(7r)R,/ T (vr ). 



(3-5b) 
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3.1. Retarded and advanced Green's functions and self-energies 

First, we consider the retarded (advanced) component of the Dyson equation (|3-1|) . 



which reads 



G 



R(A) 
E 



(vr)=G? (A) (vr) 



+ 



(d^G^- 1 (ir))G* {A) (vr) (V^G^' 1 (*)) 
vr))Gj (A) (^)(a w oG ^ ) - 1 (^)) 



(V.G^ M < 



G, 



f (i) w- 

(3-6) 



fl(ii)-i. 



W)xG„ fl(i) W(v f G; 



7T)) 



G? (A) (vr) 



(3-7) 

because of the relations given by Eqs. (|B-3j) and ()B-4|) . 

The self-energies due to the impurity scattetring are obtained from Eqs. (|2-29|) 
and flH3H) as 

d-K 1 . 



E£ '(7r u ,7r) = n imp 

AH(A)/ \ 

I 71 " i 71 ") = "-imp 



(27Tft) rf 



fiR(A), o /\ ~ 



UTt-Tt'G^ '{jT ,7T JUtt'-, 



(3-8) 
(3-9) 



in the Born approximation for the general potential, Eq. (|2-30|) . and as 

\R{A) , (K , On^t-R(^) /On 



^ {A) (VT°) 



^ {A) (vr°) 



^) / On / On^(A) /On 



(3-10) 
(3-11) 



n imp T ^ (tt u )^ (^ U )T ^ (* u ), 
in the T-matrix approximation for the <5- function potential, Eq. (|2-31|) . 

3.2. Lesser Green's functions 

Using the relation given by Eq. I)B-6|) . the lesser component of Eqs. (|3-2|) and 
(|3-3|) can be expressed as 

G^(vr) = (Gi(n) - G|(tt)) / t (tt°) 

+ Gq(tv) (S< - - £!(7r))/ T (7T )) G^vr) 



-G^vr) ((^(tt) - Z«{k))G^)(V^- 1 {k)) 



(V.Gr i (vr))G «(vr)(^(^) - JC ») ) G6>Ro/ T ( 



ft/ 



(3-12) 



and 



+ G«(vr) (T< - (27g(7r) - ^(tt))/^)) G£(tt). 



(3-13) 
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These sets of equations have the solutions of the form 

G%(tt) = G< ^ttRo/^tt ) + G< jH (k)M<k% (344a) 

= + ^.i/OO/tCA (3-14b) 

and 

<%(t) = ^(ttRo/^tt ) + G< J7 (tt)/ t (A (345a) 

^b(t) = £b>Ro/t(^ ) + JC< // (^)/ T (vr°), (345b) 

respectively. The first and the second terms represent the classical contributions from 
the states at the chemical potential with a factor of 3 7r o/^(-7r°) and the quantum 
contributions from all the occupied states with f=(n ), respectively. Substituting 
Eqs. (K-S-14al) and (I344bl) into Eq. (l3~T2l and Eqs. (I3-If>al) and (I345bj) into Eq. (l3~T3l . 
we finally obtain 

G%M = G*(tt)£< ^Gtf (tt) 



- - G«^))(V.(iY (7r) + Efa)))G$(*) 

G*(vr) (^(tt) - iV w (ffo(Tr) + i(^(»r) + £ V)))) G^(vr) 



+ 2^^(70 + G£(vr)j, (346) 

G^j/OO = G|(7r) - G|(7r), (347) 
JC< /z (7r) = ££(tt) - £|(tt), (348) 

and 

(349) 
(3-20) 
(3-21) 
(3-22) 

When Gg j(ir) is numerically calculated, Eq. ()346|) should be used as it is. 
However, when it can be calculated analytically, the calculation carried out to solve 
Eq. (|346|) is facilitated by constructing the following equation from Eq. (|346|) : 



G%! 


M 


= 0, 








= 0, 




Gb,ii 




= G^(vr) 


- Gg(vr) 


^B,II 




= ^b(vt) 


- ^b(vt) 



G% tI (ir),Ho^)\ _ + G< j(ir)^) - (vr)G< ,(tt) 



+ I ((^(vr) - r «(vr))(V.G A (vr)) - (V K (H fr) + ^(vr)))(G A (vr) - G^(vr)) 

+ (V.G«(vr))(^ A (vr) - ^(vr)) - (G a (tt) - G R (vr))(V.( J ff (7r) + £ a (tt)^ 

(3-23) 
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This is the so-called quantum Boltzmann equation for the classical contribu- 
tion^ J^ 1 ^ in the context of the Keldysh formalism or the Baym-Kadanoff equa- 
linearized in the electric field. This self-consistent equation (|3-23|) can 
also be derived by substituting Eqs. (|2-38j) and (|2-39j) into Eq. (|2-37[) . linearizing 
it in F^ 11 , substituting Eqs. (|3-14|) . and taking the classical terms proportional to 
cW T (^°). 

The classical contribution to the ^-linear lesser self-energy is given by 

j—^u^^G^jin ,™')^'--* (3-24) 

in the Born approximation for the general impurity potential, Eq. (|2-3U|) . and by 

= n^f^Yg^)^^) (3-25) 

in the T-matrix approximation for the 5- function potential, Eq. (|2-31l) . 

In conclusion, the problem at hand of obtaining the quantum distribution func- 
tions Ge(tt) and Gg(ir) is reduced to solving the self-consistent equation (|3-23|) for 
the classical contribution G^, j(tt) and the self-consistent equations (|3-6|) and ()3-7|) for 

the quantum contributions G^ A \-k) and G^ A \tt) over the occupied states, together 
with the diagrammatic representation of the self-energy according to the Feynmann 
ruleP 

§4. Second-order response 

To go beyond the linear response response to the external fields, we proceed 
to the second-order terms in Eq. (|2-3(i|) . In particular, substituting Eqs. (|2-38|) and 
()2-39j) into Eq. l)2-36j) and comparing the second-order terms in F, we obtain 



+ \ {((^W)(^4W) - (^Go 1 (vr))(a^G MV ,(vr)) 

- (fj, v) + ((fj,u) <-» (//V)) - {{nvn'v') (i/W))} 
+ \ ((^5 7rfI /G 1 (7r))(^5^/G (7r)) 

- {ji «- u) - (// <-» v') + ((w') <-» (iV)))] (4-la 

\{{(d^G^)){d^E^)) ~ (^^(tOX^GoV)) 
(/x «-» v) + {fjiv) «-► (AtV)) - (0^j"V) «-» (y'/iV))} 



+ 2 



+ \ ((^5 7rfl /G o (7r))(0^a x G o 1 (vr)) 



(/i <-> i>) — (m' <-> i/') + ((////) <-> (^')))]Go(^)- (4-lb) 
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In this section, the argument ir of the Green's functions and the self-energies is often 
suppressed. 

4.1. Magneto-transport properties 

Now we consider the magneto-transport properties as a second-order response 
expressed by the cross terms of E and B. We start from Eqs. (|4-1|) . which can be 
rewritten as 



G Ei}Bj - G 



+ ((V^) x (V.Gq) - (V ff Go *) x (V^)) .} 
+ ~ ((^oV^ 1 )) x {d^Ga) - (5^ V ff (Go ')) x (d^o V^G ) 

Q.oHEi,Bj + G-Bj^LEi + G-E^Bj 

- ((a^Go)^^.) - (d^X^oGo 1 )) 

+ ((V.G ) x (Vj^) - (V^) x (V.Go 1 )) .} 

+ i ((^oV,G ) x (^V w Go ') - (^iV^Go) x (c\oV w G ')) 



(4-2a) 



Go- 
(4-2b) 



The retarded (advanced) component can be written as 



^R(A) _ f1 



+ 



R(A) 



f,R(A) r<R(A) fiR(A)/yR(A) , ^R(A)AR(A) 
Zj E i ,B :j < ~ r +Zj Ei Bj ' Bj Ei 



\d^S^){d^G 



{((d^)(d^ A) ) - (^G^X^Gjf); 

. - (V.Gj^- 1 ) x (V„6g A >)) } 



(v.g? (a) ; 



T ((5 7r oV.G ^ ) - 1 )x(a^V 7r G5 ?(A) ) 

(a^v^Gf^- 1 )) x (s^v^gJ 



(4-3) 
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We know that G E . (tt), (jr), G b \tt) and E^ {tt) have the forms given in 
Eqs. (|3-14aj) . (|3-14b|) . (|3-15a|) and ()3-15b|) . supplemented by Eqs. (HOI, (EH), and 
l|3-16j) . Therefore, using the property given by Eq. (|B-6|) . we find that the lesser 
component of Eq. (|4-2|) also has the solution 



G%, B ^) = G%^j(n)d n of T («°) + G< ^(vr)/^ ), 



(4- 4a) 
(4-4b) 



with the self-consistent equation for G 



and U", 



G\ 



Ei,Bj,I aliKX ^Ei,Bj,Ii 
A , ^< riA 



^E^BjJ^O + ^EiJ^Bj + Zj B J (j E i ,I 

+ 1 -{{Z% - Z%){d^G$) + (Z£ - E^Gi.) 

- (d^%)(G$ - 6$) - (dAHoW + ^)){G% - G%) 

+ ((v^< tJ ) x (v^ctf) + (v w (ff (w) + x (v.g< i7 ; 

i 



i ((V.(^-i7 ( f))x(^V^ 



- V w (Fo(7r) + r ")) x (V W (G£ - G«)) 

Zj E l ,Bj,I + U Bj^Ei,I + U E l ,I 2 ^B J 

+ ~ {(Gtf - G*)^^.) + (G^. - G|)(^(^ (tt) + r A )) 

- (^g*)^. - - (a^Gf.)^ - A?) 

+ ((V^G«) x (V^< ,) + (V W G< A x (V^M + £ A )) 



(4-5a) 



~ ((V.(G 



G£))x(^V w (fr o (7r)+itf)) 



a,V,G«)x(V,(^-^)) 



^0 > 



(4-5b) 



and 



GeuBiJiW = Ge^b^) ~ G i il s i ( 7r )> 



Ei,Bj ( 7r ) - 



(4-6a) 
(4-6b) 



Equation (|4-5|) can be transformed into an analogue of the quantum Boltzmann 
equation; 



GEi Bjji ^o(tt) + G E . B jSq - £qG e . b j 
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f^< f^A fiR y-i< i y^i< Aj4 fiR f^< I f-*R / s << A< f^A 

■^E^BjJ^O ~ ^0 ^E u B h I t ^EiJ^Bj ~ ^B^E^I "+" ^B^E^I ~ ^EiJ^Bj 

+ {E A - S R ){d^G%) + (d ni G§.)(E£ - £*) 

- (G A - G R ){d^S A .) - (d^E R )(G A - G R ) 

- {G% - G%){dAH^) + A?)) " (flW^oM + ^?))(G4 " G&) 
+ ((V^< tI ) x (V^) + (V w (flbW + x (V„G< f/ ) 

(V W G *) x (V^< (J ) + (V^G< )7 ) x (V^Hoin) + A?))).} 

;V.(^ A - E R )) x (fl^Ctf) + (S^V W G«) x (V.(£ A - £ R )) 
(V^Gtf - G*)) x (a^V w (H (7r) + ^ A )) 
(^V»(ffoW + ^)) x (V W (G# - G*)) 



(4-7) 



§5. Transport properties in the single-band model 



In the single-band case, all the Green's functions G R)A,< and the self-energies 
are scalar quantities, and can be denoted as G R,A,< and S R,A ' < 1 respectively. 



Then, from Eqs. (|3-6[) and (|3-7j) . it is readily shown that Gg A (ir) = U^ a (tt) = 



-,R,A, 



and Gg A (ir) = E^ a (tt) = 0, and hence the quantum contributions from G EII to 
the electric current up to the linear response vanish. There exists only the classical 
contribution from G EI in proportion to d 7T af + (ir ), which becomes 5(ir° — fj,) at the 
zero temperature. Therefore, one has only to obtain G EI (n,p F ) and U E j(fi,p F ). 
Equation (|3-16|) or (|3-23|) is reduced to 



-,R,At 



GeAtt) 



lmG R (n) 
ImZ^(7r) L 



ZeiW ~ »V w (flb(w) +Re^(7r))l + iReV 7T G R (Tr). (5-1 



This is the Bethe-Salpeter equation for the vertex correction S EI (jr). Similarly, for 
the magneto-transport, from Eq. (|4-5p we obtain 



lmG R (7r)E<^ tI i 



7T 



i>* s ' Imi7^(7r) 
Cu ' j '(fc27< >))(0,ReC#(7r)) + (d e (H (n) + Re^(vr)))(^G< ») 



-ididdelmU^Tr^dedrnG^Tr))) 



(5-2) 



In the case of the magneto-transport, however, the quantum contributions character- 
ized by G F ' A B .(ir) and E R,A B . (ir) do not necessarily vanish unlike G e A (tt), E e ,A (tt), 

G^' (tt), and E% A (vr). From Eq. M . we can obtain the self-consistent equation 
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to calculate Gg' A B .(ir) and U B ' A B (ir), 



- ^V,(i?o( ff ) + ZS {A \*))) x (c),oV4 W W)), (5-3) 



In the following, we use the notation (e,p) instead of (ir°,7v). 

Here, we employ the self-consistent Born approximation given by Eq. (|2-40|) . 
More explicitly, we have 

dp' 



^0 (Mj P) = n imp 



(2irh) c 
dp' 



S E^bS^P) = n i 



nn " / u-hyi 

dp' 

(2irh) d 
dp' 



u(p 


-p') 


\ 2 G^,p'), 


(5-4a) 


u(p 


-p') 


| 2 G<>,p'), 


(5-4b) 


u(p 


-p') 


\ 2 G%, Bj ,iM, 


(5-4c) 


u(p 


-p') 


\ 2 Gl A Bj {^p'). 


(5-4d) 



imp J {2<Khy 

5.1. Longitudinal conductivity 

First, we calculate the longitudinal conductivity. Substituting Eq. IJ5-ljl into 
Eq. (|5-4bj) and using Eq. (|5-4d|) . we obtain 

S^ I ( f i,p)-iV p ReS^(n,p) 

= Ij^^ 2 ^m^) ^ M - lVAH " (p,) + Rer ° R(ey)) l 

= -i(r tr /T-l)V p H (p), (5-5) 
and thus 

G%j^, P ) = ~i l m( 3^ P \ -V p H (p) + zV p ReG^,p), (5-6) 

where we have assumed an isotropic dispersion Ho(p) = p 2 /2m for simplicity, and 
introduced the electron lifetime t and the transport lifetime r tr via the relations 

h/r tx = 2 J^E_\ UpF _ plF \ 2 (l- C o S e)(-lmG§( f i,p' F )) = j dDf{Q){l - cos 6), 

(5-7) 

h,T = 2 J ^jdK^A-^oi^PF)) = J dQf{Q), (5-8) 



and the differential cross-section 



f{Q) = 27rn imp m * { ^ d \u(p F -p' F )\ 2 (5-9) 
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with the polar coordinates I? = (6, ip) of p F with respect to p' F and Qd = r(d/2) /it d l 2 . 
Here, m* is the effective mass given by 



l/m* = z/mb = z(l/m + \7 2 ReU^(fi,p))\ p ^ PF 
with the renormalization factor defined by 

l/z = l-d e ReE^(e,p F )\ e ^. 



(5-10) 



(5-11) 



From Eqs. (|5-1|) and (|5-5|) . the conventional expression for the longitudinal elec- 
tric conductivity at the zero temperature is reproduced as 



2m 
e 2 H 



dp 



r{d iH Q (p))G% U,p)d E f + {e) 



dp 



{2Tth) dK '" 

(5-12 



n*e 2 T tI m* n*e 2 r^ 
Oij = Ojj o 
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in the absence of the magnetic field. We have introduced n* = zn with the electron 
density 

n = (2 d p d F /d{2irh) d (5-13) 



and the quasiparticle transport lifetime 



'tr 



Ttr/Z. 



(5-14) 



It is important that the renormalization of the velocity due to the self-energy disap- 
pears in the above expression for an, in agreement with the Kubo formula. 

5.2. Hall effect 

We proceed to the second-order response of the electric current to the cross term 
of the electric field and the magnetic field. In the present model, using Eqs. (|5-5J) 
and (|5-6j) . Eq. (|5-2j) is reduced to 

G E t ,B 3 ,i( £ iP) ~ itijej^{d p eReG${e,p)) 

= G$(e,p)G$(e,p) [^. B . j7 (e,p) - ie ije ^d pe (H (p) + ReZ$(e,p))] , (5-15) 
as shown in Appendix O Therefore, from Eqs. (|5-4a|) and (|5-4c|) . we obtain 
ZeuB^A^P) ~ ieiji^d p eReU^(fi,p) 



dp' 



Ttr 



-ieije(T tT /T - l) — (d pt H Q (p)) 



(5-16) 
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and thus 



lmG§(e,p) r tr 



(dpeHoip)) - (fi^ReG^p)) 



(5-17) 



_Imi7Q ? (e,p) r 

as also shown in Appendix [0 Consequently, we obtain the Fermi-surface contribu- 
tion to the Hall conductivity as 

.3^2 



6 /* /* dp 

a v =Bi ^i j d£ J j2^^ dpiHo ^ G ^> B ^ £,pW ^ 



e 3 h 2 



2iri 
e 3 H 2 



dp 

J^jd (p) )G< , iBtJ (/*, p) 

dp r tr 



2vr 7 (2irh) d hm 



x (d pi H (p)) Gq (n,p) Gq (fj,, p) — (d p iHo(p)) - (^ReG^p)) 



£ijeB 



* 3 2 

n e Ttj. m. 



D £ 2 1 ^ D ^ 2 ^ 

C^-D (T.-.- — : — €inPt> 0*;, 



, , (5-18) 

Next, we proceed to the quantum contribution determined by G B ' A B (e, p). From 
Eq. (|5-3|) . we obtain 

Gg&fop) = G? A) %.P)^,(e,P) 

" X ((^( ff o(P) + ^ {A \e,p))){d £ d pl G^ A \e,p)) 



+ (^^(tt))(#C^(*)) 



2 ^(A), 



(5-19) 



In particular, in the case of a ^-functional potential, dpiS^ vanishes, and so it is 
reduced to 

(5-20) 



^ A \e,p)d e d pl G^ A \e,p) 



" ' Am 

Thus, we obtain the quantum contribution to the Hall conductivity as 

e 3 h 2 



a 



11 



Bp 



-B 



2m 
e 3 h 2 



de 



dp 



de 



B e eijp 



2tt 

e 3 h 2 1 



. )7Th)1 (d p ,H (p))(-2iImG^ :Be (e,p)f + (e)) 
^^(^Fo(p))Im^( e ,p) 9£ / + ( £ ) 



dp 



lmG*(n,p). 



(5-21) 



2vr 2m 2 J (2nh) d 

This term is negiligible in the usual metal compared with the classical contribution 
given by Eq. (|5T8|) . The non-local self-energy produces a small correction to ths 
expression, which is of higher order in 1/r. Therefore, it is negligible in the weak 
scattering limit. 

From the above results, we reproduce the conventional expression for the Hall 
coefficient 



R = Pyx /B 2 



1 m 



n e m 



1 

neb 



(5-22) 
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The above calculation gives an expansion of the Hall conductivity /resistivity 
in co c T tT <C 1 with the cyclotron frequency lo c = eB/m and a weak magnetic field 
B. We also note that the above results can be easily extended to the case of a 
finite temperature T by exploiting the Sommerfeld expansion and expanding the 
self-energies and the vertex corrections in T to obtain the coefficients to the T 2 
terms in the conductivity tensor and the Hall coefficient. 

§6. Discussion and conclusions 

We have developed a new theoretical framework to study the electromagnetic re- 
sponse of the charged system in terms of the Keldysh formalism. The Moyal product 
and its expansion are the key concepts in this approach, and the gauge invariance is 
implemented from the beginning of the calculation. Here we discuss the relationship 
of the present approach to others, which are equivalent, in principle, within their 
domains of validity, but differ with regard to procedures of their calculations. 

The most well-established formalism is the Kubo formula for the linear response, 
where the retarded response functions are obtained from the Matsubara Green's 
function via the analytic continuation, which can be calculated using the diagram 
method. In this formalism, it is necessary to take the diagrams for the self-energy and 
vertex correction consistently in order to satisfy the Ward-Takahashi identity. This 
is rather straightforward for the conductivity, but it requires an ingenious treatment 
for the Hall conductivityP In the latter case, to guarantee the gauge invariance, the 
factor q 2 5 a b — q a q\, where q a is the wavenumber of the vector potential of the external 
magnetic field must be extracted and a careful treatment of the self-energy and ver- 
tex corrections is required, as mentioned above. In our formalism, by contrast, the 
equation of motion is gauge-invariant from the beginning, because of the replace- 
ment of the argument, p — > tt = p — qA, and the fact that the vertex correction is 
represented by G < and E < , while the self-energy by U R ' A . It is noted here that both 
E < and E RA satisfy a single equation represented in the Keldysh space. Another 
advantage of the present approach is that no analytic continuation is needed. This 
is particularly useful when one analyzes non-linear responses, which is beyond the 
scope of the Kubo formalism. 

Another standard treatment of the non-equilibrium state is given by the Boltz- 
mann transport theory. This is closely related to the Keldysh formalism in spirit, 
and G K is the quantum analogue of the distribution function f(k, r) appearing in the 
semi-classical Boltzmann equation. Actually, the quantum Boltzmann equation is de- 
rived from the Keldysh formalism, and by assuming Gq(tt) = 27ri5(-7r — e(7r))/ T (-7r ) 
with e(n) being the energy dispersion of the particle, it reduces to the semi-classical 
Boltzmann equation. This formalism is akin to the present one, but the change of the 
self-energy E R)A and the Green's function G R ' due to the external electromagnetic 
field is often ignored in the linear response regime. This semi-classical assumption 
is valid for single-band models, but not for multi-band systems, since in the latter 
case, the quantum contribution is missing in this assumption. Our formalism can 
also treat this multi-band case systematically. 

Lastly, it is interesting to note the relation of the present formalism to the Stf eda 
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formula^ for the linear response. Stfeda decomposed the electric conductivity ten- 
sor an into a/.- + erf/, where the classical contribution cr/ - comes from fermions at 
the Fermi surface while the quantum contribution aH from all the occupied states. 
Without the magnetic field, Eq. (gHH) with Eqs. (j^B|) . (EUHl), and (j3T7|) coincides 
with the Stfeda formulcP^ in the absence of tj_. Even under a magnetic field, a similar 
decomposition of Eq. (|4-4j) facilitates the calculation of the conductivity, as explicitly 
shown in ^15.21 This decomposition is particularly useful when one investigates the 
intrinsic and extrinsic contributions to the Hall conductivity, an interpretation of 
which has been controversial for anomalous Hall effect for many decades but clari- 
fied in a unified wayP" Our present formalism provides a method to calculate err- 
and ajj diagrammatically even in the presence of the particle-particle interaction. 
Therefore, the present approach can be regarded as a further development of the 
Stfeda formula. 
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Appendix A 

Moyal Product 



We start from the Wigner mixed representation of T and £_ in terms of (X;p): 
t(X;p) = J dxe^^f (^±£, , (A-la) 

mX;p) = J dxe^^f (*±£, . (A-lb) 

Then, the Fourier transform of the convolution given by Eq. ()2-lUj) reads 
(I * lm P)=jdx j dx 3 f (^±£, § + **) t (f + *3, ^) 



dp' jjjsr-.vsvn f dp" ( _„ v ,., :! __ r . 



dx d x, / «P e *M#-«S)/» 



(2iTh) d + 1 ' J (2irh) d + 1 

dX > [ dX" ( T ^-- f T ^— T e V'^v'^),n 

J J (2-Kh) d + l J (2TTh) d+l 

(X' \ ( X" 

x± \X + —;p + P '\S\X-—; P - p" 

dX > [ dX" I -^-r I eW^p'^yn 



(2irh) d+1 J (27rh) d+1 

dX' f dX" 



dp 1 f dp" 



(2lTh) d + 1 J (27Th) d + 1 
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= f(X-,p)e^ x ^^-* d ~^ x ^E{X-,p). (A-2) 
Appendix B 

Basic Properties of the Retarded, Advanced, and Lesser Components of 

Product of the Matrix-Form Functions in the Keldysh Space 

Given the matrix-form functions in the Keldysh space, 

with i = 1, • • ■ , n, their product 

t(p)=ti(p)t 2 (p)---tn(p) (B-2) 
satisfies the following relations: 

t R (p)=f R (p)f R (p)...f R (p), (B-3) 
f A (p) = f A {p)f A (p)---f A {p). (B-4) 

Furthermore, when each T_ i satisfies the equilibirum condition 

f<(p) = (f A (p)-f R (p))f T (e), (B-5) 
then the lesser component of the product is reduced to 

f<(p)=f<f A ...f A 

+ fRf<fA ...fA 

| 'yRryRryK 'y-A sy^A 

+ ■•• 

! y-R y-R -y , < y~~A 

+ fR...fR_ i f< 
= {f A (p)-f R (p))f T (e) (B-6) 

for n > 2. 

Appendix C 

Derivation of Eq. \5-17\) 

We start from Eq. (|5-2[) . which has been reduced from Eq. ()4-7|) in the single- 
component scalar problem. Then, substituting Eqs. (|5-6|) and (|5-5|) into Eq. 1)5-2(1 
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leads to 



hnG$(e,p)E<Ae t p) 



+ i 



'1 {(d/-±d pi H Q {p)){d^G*{e,p)) 
+ (d pe (H (p) + ReU^e,p)))(d pel ^G§(e, P )G^(e,p)d pl H (p))} 



"f- {(d pl d pl (H (p) + ReZK(e,p)))(d pe ,ReGK(e,p)) 
+ (d pe (H (p)+ReZK(e,p)))(d pe ,d p >ReGK(e,p)) 
- d pi ((a p ami7 ii (M ! p))(^ImG^( At ,p)))} 



ImGX(E,p)Z<r(e,p) 



ImE*(e,p) 
+*^{(d pe ^d p >H (p))(d pe/ ReGK(e,p)) 

- (d pe ^d pl H (p))G§(e,p)Gi(e,p)(d 7Te/ (H (p) + ReE^(e,p)))) 



.T tr 



(lmZK(e,p))(d p >H (p))(d p tG«(e,p))(d,GS(e,p)) 



Hi 



d pl Re [(d pi (H (p) + S^,p)))(d pel G^,p))) } 



(C-l) 



The last line of the above equation identically vanishes due to the anti-symmetric 
tensor ewj. Using eu'j{d p igi{p)){d p e< 52 (p)) = 0, Eq. (|C-1|) can then be rewritten as 
Eq. (|5-15l) in the present model Hq(p) = p 2 /2m. 

Multiplying both sides of Eq. (|5-15|) by |«p- p '| 2 ^/ r tr and performing a momen- 
tum integration, we obtain 



II r 



dp' 

dp' 



(2-kK) 



d \ u p-p / 



^G%^p') + l e -^{d p>l ReG^p')) 
. 2 1mG$(e,pf) 



ImZK(e,p>) 



^ Bi Ae,p') + i e -^(dMH (p')+Re^(e,p'))) 



m 



(C-2) 



Using Eqs. (|5-4a|) and (|5-4c|) . we reproduce Eq. (|5-16|) . 
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